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ABSTRACT 



Context. After the detection of a 321-day periodicity in X-rays, HESS J0632+057 can be robustly considered a new member of the 
selected group of gamma-ray binaries. These sources are known to show extended radio structure on scales of milliarcseconds (mas). 
Aims. We present the expected extended radio emission on mas scales from HESS J0632+057. 

Methods. We observed HESS J0632+057 with the European VLBI Network (EVN) at 1.6 GHz in two epochs: during the 
January/February 201 1 X-ray outburst and 30 days later 

Results. The VLBI image obtained during the outburst shows a compact ~0.4 mJy radio source, whereas 30 days later the source has 
faded and appears extended, with a projected size of ~75 AU. The peak of the emission is displaced between runs 21 ±5 AU, which is 
bigger than the orbit size. The position of the radio source is compatible with the Be star MWC 148, which sets the proper motion of 
the binary system below 3 mas yr"' in each coordinate. The brightness temperature of the source is above 2 x 10* K. We compare the 
multiwavelength properties of HESS J0632+057 with those of the previously known gamma-ray binaries. 

Conclusions. HESS J0632+057 displays extended and variable non-thermal radio emission. Its morphology, size, and displacement 
on AU scales are similar to those found in the other gamma-ray binaries, PSR B1259-63, LS 5039, and LS 1 +61 303, supporting a 
similar nature for HESS J0632+057. 
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1. Introduction 

Gamma-ray binaries are binary systems containing a compact 
object orbiting a companion star, which are able to produce 
GeV and/or TeV emission and show the peak of the spectral en- 
ergy distribution (SED) at MeV-GeV energies. The three "clas- 
sical" gamma-ray binaries are PSR B1259-63, LS 5039, and 
LS I -1-61 303, which have been unambiguously det ected as 
point - like sources at TeV energies (s ee the reviews in lParede"sl 
I2008t iBosch-Ramon & Barkovll201 lb . The broadband emission 
from radio to very high energy (VHE) gamma rays of these 
sources is variable and periodic, following the orbital period of 
the binary system. These three sources present variable milliarc- 
second (mas) scale radio structures that have been detected with 
Very Long Baseline Interferometry (VLBI). We note that other 
binary systems have been detect ed up to GeV energies, but not 
confir med at TeV energies (e.g., iTavani et alJl2009 : lAbdo et~an 
I2OO9I for Cy gnus X-3), inclu ding the garnma-ra y binary can- 
didate IFGL J1018.6-5856 (ICorbet et aUlIoTlt iPavIov et alJ 
I2OII '). Two main scenarios have been proposed to understand 
the multiwavelength behaviour of gamma-ray binaries, the basic 
difference being the nature of the particle accelerator. In one of 
them, particles are accelerated in the jets of a microquasar (see, 
e.g., Bosch-Ramon & Khangulvan 2009). In the other one, parti- 
cles are accelerated in the shock between the relativistic wind of 
a young non-accreti ng pulsar and the stellar wind of the massive 
companion star (see lTavani & AronllI997l : lDubusll2006l) . 

HESS J0632+057 is the most recent addition to the selected 
group of gamma-ray binaries. The HESS Collaboration found 



the p oint-like TeV source HESS J0632-H0 57 dAharonian et al.l 
2007), which was variable at TeV energies dAcciari e t al. 200^ 
The s ource has a varia ble X-ray counterpart dHinton et aU 
I2OO9I: lAcciari et alj II2009I) . as well as a variable radio counter- 
part (ISkilton et al.ll2009h . The massive BOpe star MWC 148 
(HD 259440), located at ~1.5 kp c, was proposed as the opti- 
cal counterpart dHinton et alj|2009l) . The SED from radio to TeV 
energies is very simila r to the one of LS I -1-6I 303, which also 
contains a Be star ( Hin ton et a n i2009l) . All these results suggest 
that HESS J0632+057 is a new gamma-ray binary, displaying an 
SED one order of magnitude fainter than LS I +61 303. This is 
an important fact that could shed light on the luminosity distri- 
bution of this new population of binary systems in the Galaxy. 

Recently, a periodicity of 321 +5 d has been re- 
vealed thanks to long-term X-ray observations conducted with 
SwiftPCRT, str ongly supporting the binary nature of the source 
dBongiorno et al. 2011). This orbital period is between the pe- 
riod of PSR B1259-63 (1237 days) and the one of LS I +61 303 
(26.5 days). The X-ray light curve, which covers three orbital 
cycles, shows a bright peak and a fainter secondary peak sep- 
arated ~0.5 in phase. On the other hand, optical spectroscopic 
observations have not allowed the binary nature of the source t o 
be unveiled up to now (Arag ona et aIj2010l : ICasares et al l20Tlh . 

The last bright X-ra y outburs t of HE SS J0632+057 occurred 
in February 2011. Falc one et al.l d201 Ih report increased X-ray 
activity detected by Swift/XRT between January 23 and February 
6, 2011 (MJD 55584-55598). The VERITAS Collaboration re- 
ports increased activity at energies above 300 GeV betwee n 
February 7 and 8, 2011 (MJD 55599-55600) dOnil20Ilh . 
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The MAGIC Collaboration reported increased gamma-ray flux 
above 200 GeV during February 7-9, 201 1 (MJD 55599-55601) 
(Mariotti 2011), confirming the VHE active state and lowering 
the measured energy threshold. During the outburst, we observed 
the source with the European VLBI Network (EVN) to explore 
its radio emission on mas scales. 

The non-thermal radio emission produced by gamma -ray bi- 
naries has been observed on scales of 1-100 mas. M assi et alj 
( |2004|) report a jet-like and precessing structure with an exten- 
sion of 10-80 mas (25-200 AU) from LS I +6 1 303. VLBA 
obser vations at 8 GHz along the orbital cycle (iDhawan et alj 
12006 ) have shown that the ~7 mas (14 AU) structure follows 
a cometary behaviour, while the peak of the emission traces 
approximately an ellipse. For LS 5039, extended radio emis- 
sion has been found at 5 GHz, with projected sizes of 6 mas 
(15 AU) with the VLBA (Paredes et al. 2000; Ribo et al. 2008), 
and betwe en 60 and 300 mas (150-750 AU) with the EVN and 
MERLIN dParedes et all 12002) . Observations with the VLBA 
covering an orbital cycle show periodi c morphological vari - 
ability and hints of peak displacement C Moldon et alJl20l"T3) . 
Recently, Moldon et al. (201 la) has shown fliat PSR B1259-63, 
the only gamma-ray binary that is known to contain a pulsar, dis- 
plays emission with a projected extension of 50 mas (120 AU) 
at 2.3 GHz during periastron. The peak of the radio nebula is 
detected at distances between 10 and 50 AU from the system. 

Here we present the first VLBI radio images of 
HESS J0632+057, obtained with the EVN during and af- 
ter the high-energy activi ty period in early 20 11 (preliminary 
results were presented in iMoldon et alj [201 Ibh . The observa- 
tions and data reduction are described in Sect.|2] The results can 
be found in Sect. [3] A discussion of these results in the context 
of gamma-ray binaries can be found in Sect. H] 

2. Observations and data reduction 

Following the report of X-ray activity fro m HESS J0632-F05 7 
between January 23 and February 6, 201 1 (iFalcone et alJl20lll) . 
we observed the source with the EVN in target of opportunity 
(ToO) mode in two epochs separated by 30 days. The first ra- 
dio continuum observation (run A) was conducted at 1.6 GHz 
(18 cm) during eight hours on February 15, 201 1 (MJD 55607), 
UTC 15:50 to 00:00. The antenna network consisted in seven 
antennas: Effelsberg, Jodrell Bank Lovell, Medicina, Onsala, 
Torun, Westerbork, and Hartebeesthoek, providing baselines 
ranging from 200 to 8 000 km (200 to 1000 km without 
Hartebeesthoek). A data rate of 1024 Mbps per station was 
directly streamed to the central processor at JIVE and corre- 
lated in real time (e-VLBI). The rapid response of the e-EVN 
observations and correlation allowed us to perform a second 
ToO (run B) in a disk-recorded session with 12 EVN anten- 
nas: Effelsberg, Jodrell Bank Lovell, Medicina, Onsala, Torun, 
Westerbork, Hartebeesthoek, Svetloe, Badary, Zelenchukskaya, 
Nanshan (Urumqi), and Sheshan (Shanghai), providing base- 
lines ranging from 200 to 10000 km. The observation lasted 
ten hours and was conducted on March 17, 201 1 (MJD 55637), 
UTC 13:00 to 23:00, with a data rate of 1024 Mbps. This second 
epoch provides considerably better sensitivity. The EVN project 
codes for the observations are RR005 and RM006, respectively. 

Both observations have a similar structure, switching be- 
tween HESS J0632-H057 and the phase reference calibrator 
J0619-1-0736, separated 3^9, with cycling times of five minutes 
to avoid losing phase coherence. The phase reference calibra- 
tor was correlated at the position Q'j200().() - 06^19™09';9710 and 
(5j2()OO.o = 07°36'41'.'220 in the frame of ICRF (Goddard VLBI 



global solution 2010^. This position has absolute uncertainties 
of 2.0 and 1.1 mas in right ascension and declination, respec- 
tively. The fringe finders DA193, 3C147, and OQ208 were ob- 
served during run A and 0528-1-134, and DA 193 during run B. 

The data reduction was performed in AIPSQ. Flagging based 
on predicted off'-source times, owing to slewing or failures, was 
applied using UVFLG. A priori visibility amplitude calibration 
used the antenna gains and the system temperatures measured at 
each station. The amplitude calibration was improved by scal- 
ing the individual antenna gains by a factor obtained from the 
phase calibrator and fringe-finder models. We used ionospheric 
total electron content (TEC) models based on GPS data obtained 
from the CDDIS data archiv^ to correct the global phase vari- 
ations due to the ionosphere. The parallactic angle correction 
was applied with VLBAPANG. Fringe fitting on the phase refer- 
ence calibrator was performed with the AIPS task FRING, and 
the solutions were applied to the target source. The instrumental 
bandpasses were corrected using BPASS. The data were aver- 
aged in frequency and time, and clean images were produced 
with IMAGR. A cell size of 1 mas was used for cleaning both 
epochs. The images were produced using a weighting scheme 
with robust parameter for run A, and 2 (which is slightly more 
sensitive to extended emission) for run B. A tapering of 30 MA 
was applied to avoid the presence of possible unreliable high- 
resolution features due to sidelobes of the synthesized beam. No 
self-calibration of the data was possible because of the low flux 
density of the target source. 



3. Results 

The resulting VLBI images at 1 .6 GHz are shown in Fig. [T] 
Run A shows a 410 + 90 yuJy point-like source. Run B, observed 
30 days later, shows a 180 + 30 /iJy (measured with TVSTAT 
within AIPS) faint source with extended emission towards north. 
The Icr rms noise close to the source is 50 /iJy beam"' and 
13 fi}y beam"', respectively. We note that the flux density of 
the extended emission detected in run B is below the rms noise 
limit of the image of run A. In Table [T] we show some parame- 
ters of the observations, the obtained synthesized beam, and the 
fitted components to the source (Gaussian components obtained 
with JMFIT within AIPS). 

The emission in run A is described by a component located 
at Q-jiooo.o = 06''32™59^.2567(1) and ^jjooo.o = 05°48'01'.' 162(2). 
The values in parenthesis refer to the uncertainty in the last digit 
(1.3 and 2 mas, respectively), and correspond to the formal er- 
rors of the fit. Run B is described well by a point-like core (see 
Table [U and an extended component (which has an intrinsic PA 
of 4 + 20°) located at 21+6 mas in PA 30 + 10° with respect to 
Core2. The total extension of the source is ~50 mas. 

The relative astrometry shows that the peak of the emission 
suffers a total displacement of 14±3 mas in PA 190 + 9° between 
runs A and B. This displacement is aligned with the direction 
of the extended emission. The separation between the peak po- 
sitions corresponds to a projected linear distance of 21+5 AU. 
Covering this distance in the 30 days between runs requires a 
constant velocity of 1200 ± 300 km s"'. 



' http://gemini.gsfc.nasa.gov/solutions/2010a/2010a.html 

- The NRAO Astronomical Image Processing System. 

http://www.aips.nrao.edu/ 
^ The Crustal Dynamics Data Information System 

http://cddis.nasa.gov/ 
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Fig. 1. EVN radio images of HESS J0632+057 at 1.6 GHz during the X-ray outburst and 30 days later North is up and east to the 
left. The observation dates are quoted at the top of each panel. The synthesized beam (see Table [TJ is displayed in the bottom-left 
corner of each image. The red crosses mark the position and 3cr uncertainty of the fitted component of run A. For each image, 
the displayed contours start at 3cr and increase by factors of 2'''^; the Icr rms close to the source is 50 yuJy beam"' in run A and 
13 /iJy beam ' in run B. 



Table 1. Observation parameters and Gaussian components fitted to the images. 



Run 


MJD 


Num. 


^HPBW 




Comp. 


^ peak 


S 1.6 GHz 






Maj. Axis 


Min. Axis 


Type 






anten. 


[mas^] 


[°] 




[j"Jy b-'] 




[mas] 


[mas] 


[mas] 


[mas] 




A 


55607.83 


7 


36x23 


2 


Corel 


340 ± 50 


410 ±90 






12 ± 13 


12 ± 10 


Point-like 


B 


55637.75 


12 


35x19 


-7 


Core2 


81 ± 14 


90 ±25 


-3 ±2 


-14 ± 3 


9 ± 11 


12 ± 12 


Point-like 












Extended 


56 ± 13 


110 ±40 


6±3 


5 ±6 


48 ± 18 


11 ± 10 


Extended 



"Displacements with respect to the red crosses in Fig.[T] 



4. Discussion 

The first VLBl image of HESS J0632-I-057 at 1.6 GHz, taken 
a few days after the peak of the January 2011 X-ray outburst, 
shows a compact 410 juJy point-like source (see Fig. [1]). After 
30 days, the source displays one-sided extended emission with 
a total flux density of 180 yuJy and with an extension of 50 mas 
(75 AU assuming a distance to the system of 1.5 kpc) in PA 30°. 
The peak of the emission is displaced 14 mas (21 AU) in PA 
190°. We note that the rms noise of the first image is above the 
flux density of the extended emission in the second image. 

Even though the orbital parameters of HESS J0632-I-057 are 
unknown, the size of the orbit can be derived from the orbital 
period (321 days), assuming a mass of the co mpact object of 
1.4 Mo, and a Be star of 16+3 Mq ( Aragona et al. 2010). The 
semimajor axis of the orbit is 2.4+0.1 AU, which corresponds 
to a projected distance of 1.6+0.1 mas (or less considering the 
inclination and the argument of the periastron of the orbit). The 
peak displacement reported here, 14 mas, cannot be explained 
by any motion within the orbit. 

The measured position in run A is clearly compatible with 
the UCAC3 catalogue position of the Be star MWC 148 , which 
has an uncertainty of 14 mas in each coordinate (Zacharia s et alj 
12010 ). therefore the detected radio source is unambiguously re- 
lated to the Be star We computed the proper motion of the source 
using the position from run A with an uncertainty of 25 mas to 
include possible peak displacements, the position of M WC 148 
from UCAC3, obtained at mean epoch 1984.2 (Zachari as et al. 
120101). and the VLA-D and GMRT positions from Skilt on et al. 
(l2009l) (with uncertainties of 300 and 500 mas, respectively). 



Despite the time span of 27 years of observations, the proper 
motion is compatible with no motion. The result is fia cos{6) = 
0.2 ±1.0 mas yr"' and fig - -0.2 ± 0.9 mas yr"'. 

The brightness temperature corresponding to the radio 
source in run A (see Table[T]l at 1.6 GHz is 2x10* K. Such a high 
value, together with non-simultaneous radio observations sug- 
gesting a negative spectral index around ~ -0.6 ( Skilt on et al] 
2009), rules out a thermal emission mechanism. Non-thermal 
synchrotron radiation remains, therefore, as the most plausible 
interpretation for the HESS J0632+057 radio emission. In this 
case, and assuming equipartition between the relativistic elec- 
trons and the magnetic field, the observed properties correspond 
to a total energy between 0.1 and 100 GHz of 4x10^^ erg, with 
an equipartition magnetic field of ~0.02 G. These values are ap- 
proximately ten times lower than those obtained for LS 5039 
from VLBl observations dParedes et ai.ll2000l) . 

The gamma-ray binaries PSR B 1259-63, LS 5039, and 
LS 1 +61 303 show a radio morphology with a central core 
and one-sided extended radio emission on scales of a few AU 
(although bipolar exten ded emission has a l so been detected 
at some orbi tal pha ses: iDhawan et al] 120061: iRibo et alj 120081: 
iJVIoldon etalJl2011all3) . In all these cases, the authors find mor- 
phological changes on timescales close to the orbital period, and 
displacements of the peak of the emission bigger than the orbit 
size. The morphology and peak displacement reported here show 
that the radio emission of HESS J0632+057 is similar to the ones 
observed in the other known gamma-ray binaries, suggesting a 
similar nature for HESS J0632+057. 

In Fig. |2] we show the phase-fold ed light curve of the X- 
ray data from lBongiorno et al.l (1201 ih . together with radio (top 
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Fig. 2. Multiwavelength light curve of HESS J0632+057. The 
data were folded with a period of 321 days. Zero phase has been 
arbitrarily defined at MJD 54857. The grey circles are X-ray 
observations by SwiftpCRT (0.3-10 keV) from lBongiorno et alj 
(|201 11) . Top: the red circles are the EVN observations from this 
paper. The other radio data are from Skiltonet al. (2009). The 
black circles (filled in grey) are Swift/XRT data obtained during 
the same cyc l e as t h e EVN data . Bottom: TeV emission from 
lAcciari et al.1 (l2009l) : \On^ ( 1201 lb : iMariottil (1201 ll) : iMukheried 



panel) and TeV (bottom panel) data. Our two EVN measure- 
ments suggest a radio flux decrease of the most compact re- 
gion of the source similar to the X-ray one after the outburst, 
although at a slower rate. However, a VLBI monitoring of the 
outburst is required to understand the flux changes better On the 
other hand, a persistent and higher flux was measured in previous 
lower resolution GMRT observations at a similar frequency up to 
~100 days after the outburst. This can be explained by the pres- 
ence of diffiise (resolved) emission, which is lost on VLBI scales 
(the shortest significant baseline in our data samples spacial fre- 
quencies around 1 MA, which corresponds to angular scales of 
200 mas). 

LS I +61 303 shows an X-ray out burst before apastron 
that is correlated wi th the TeV outburst (lAnderhub et al.ll2009t 
IZabalza et al.l 1201 ll) . The duration of the X-ray outburst of 
HESS J0632H-057 is ~5 times longer (although ~0.1 in phase in 
both cases). We computed non-simultaneous TeV/X-ray flux ra- 
tios during the 201 1 outburst (around phase 0.31) and compared 
them with the ones of LS I -h61 303. The TeV/X-ray flux ratio 
of HESS J0632+057 is ~5 higher (the X-ray flux is lower, or the 
TeV emission is higher) than the one of LS I +61 303. This is 
consistent with a lower magnetic field in HESS J0632-I-057. On 
the other hand, TeV emission has been detected at similar levels 
during the X-ray dip, and thus does not follow the X-ray light 
curve (Fig.|2|. When comparing the average radio (0.7 mJy) to 
X-ray flux ratios, we see that the X-ray flux is higher, or the 
radio emission is lower, in HESS J0632-I-057 by a factor ~30. 
The outbursts from these two sources are observationally very 



different, although the lack of information regarding orbital pa- 
rameters prevents any further discussion. 

As a final comment, we could interpret the displacement be- 
tween the peaks of runs A and B (14 mas in 30 days) as produced 
by the proper motion of a blob ejected within a microquasar (see 
Mirabel & Rodriguez 1994forGRS 1915+105). A blob velocity 
above 0. Ic requires that the jet is pointing towards the observer 
with an angle below 2°, which would be either very restrictive 
or fine tuned. 

In conclusion, HESS J0632+057 displays extended and vari- 
able radio emission at 50-100 AU scales, with a projected dis- 
placement of the peak of the emission of 21 AU in 30 days. 
Similar morphologies and behaviours have been found in the 
other gamma-ray binaries. However, a more detailed monitoring 
of the variability of the source with high-sensitivity VLBI ob- 
servations along the orbital cycle is required to measure the total 
extension and morphology of the extended emission, to mea- 
sure the morphologic and astrometric changes at different orbital 
phases, and to search for periodicity in these changes. 
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